UFTP preprint 407/1996, TUM/T39-96-1 



IR-Renormalon Contribution to the 
Longitudinal Structure Function Fl 

E. Stein^, M. Meyer-Hermann''', L. Mankiewicz''^ , and A. Schafer^ 

^Institut fiir Theoretische Physik, J. W. Goethe Universitat Frankfurt, 
Postfach 11 19 32, D-60054 Frankfurt am Main, Germany 

^^Institut fiir Theoretische Physik, TU-Miinchen, D-85747 Garching, Germany^ 



o 

cn 

(N 
> 

in 

m 
T— I 

o 

^ , Abstract: The available data on Fl suggest the existence of unexpected large higher twist con- 
Ph! tributions. We use the 1/A^/ expansion to analyze the renormalon contribution to the coefficient 
o I function of the longitudinal structure function -Ff ". The renormalon ambiguity is calculated 
^ ■ for all moments of the structure function thus allowing to estimate the contribution of "genuine" 
^ ■ twist-4 corrections as a function of Bjorken-x. The predictions turn out to be in surprisingly 

• ^ ', good agreement with the experimental data. 
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One of the interesting quantities that can be measured in deep inelastic lepton-nucleon scat- 
tering is the ratio 

of the total cross-sections for the scattering of longitudinally respectively transversely polarized 
photons and a nucleon, where M is the nucleon mass and Fl = F2 — 2xFi. This ratio provides 
a clean test of the QCD interaction since it vanishes identically in the naive parton model. The 
experimental information on R{x, Q^) is still limited but much better data should be available 
in a few years from now. Phenomenological fits to the existing data p| suggest surprisingly large 
higher-twist corrections of the form 



in(gvA2) \^ g2 + 1 0.1252 + xV g2 Q4 + 0.3 

where A = 0.2 GeV, hi = 0.0635, 62 = 0.5747 and 63 = —0.3534 and all momenta are in GeV. So 
far the leading perturbative corrections and the target mass corrections to R have been calculated 
0. The genuine power suppressed, i.e., l/Q"^ corrections can be analyzed in the framework of 
Operator Product Expansion ^. The power suppressed corrections which still have to be 
determined arise from matrix elements of higher twist operators and are sensitive to multiparton 
correlations within the target. An estimate of these matrix elements is a delicate problem which 
could not yet been solved. 

By comparing the experimental data with the known corrections, it was possible however 
to disentangle mass corrections and true higher twist corrections ^ and even to estimate the 
magnitude of twist-4 matrix elements contributing to the second moment of the nucleon structure 
function F^ and F2 

In the present paper we shall use the one-to-one correspondence between ultraviolet renor- 
malons (UV) in the definition of higher twist corrections and infrared renormalons (IR) in the 
perturbative series which defines the twist-2 contribution to investigate the structure of power- 
suppressed corrections to the longitudinal structure function F^. 

IR-Renormalons have recently received much attention because of their potential to generate 
power-like corrections. For a physical quantity like F^ the perturbative QCD series is not 
summable, even in the Borel sense, due to the appearance of fixed sign factorial growth of its 
coefficients. It results in a power-suppressed ambiguity of the magnitude ~ ^qcd/Q^ [B- Such 
terms show the need to include higher twist (non perturbative) corrections to give a meaning 



to a summed perturbation series 0, [T^. On the other hand also the higher twist corrections 
themselves are ill-defined. The ambiguity in their definition, due to the UV-renormalon, cancels 
exactly the IR-renormalon ambiguity in the perturbative series which describes the twist-2 term. 
In turn, the investigation of the ambiguities in the definition of the perturbation series of leading 
twist shows which higher twist corrections are needed for an unambiguous definition of a physical 
quantity. 

In practice, one has observed the empirical fact that in cases where the perturbative series was 
studied in parallel with the higher twist corrections, such as the polarized Bjorken sum rule and 
the Gross Llewellyn-Smith sum rule, the ambiguities produced by IR renormalons in the leading 
twist contribution were roughly of the same order of magnitude as the best available theoretical 



estimates of the higher twist corrections |]TT|. Thus, despite fundamental objections |jT2|, for 



phenomenological purposes one may use IR renormalons as a guide for the magnitude of higher- 



twist corrections [13|. The obvious advantage of such an approach is that the IR calculation 
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can be done for all moments, and hence the result can be extended to the full x-dependence 
of the higher-twist contributions. One has to keep in mind, however, that the last step is even 
less justified, as the order of magnitude correspondence between IR ambiguities and higher twist 
corrections has been tested only for sum rules for first moments of structure functions. 
We focus on the flavor non-singlet part of the longitudinal structure function 



(3) 



i.e. on the difference between the proton and neutron structure function, and calculate the 
infrared renormalon contribution. This will also provide the exact coefficients of the perturbative 
series of Fl in the large Nf approximation |T^, [1^, |T^. In the framework of the 'Naive Non 
Abelianization' |17, 18, 19, 2^ this can be used to approximate the non-leading Nf terms. 

We start with the well known hadronic scattering tensor of unpolarized deep inelastic lepton 
nucleon scattering parameterized in terms of two structure functions F2 and F^. 



Zn J 



2x 



{p-qf 



P^,qu + Puq^, 
p-q 



(4) 



Here is the electromagnetic quark current, x = Q"^/ {2p ■ q) and q^ = —Q^. The nucleon state 
\p) has momentum p (averaging over the polarizations of the nucleon is understood). The non- 
singlet moments of the structure functions F^~^ {k = 2, L) can be expressed through operator 
product expansion M in the following form: 



N-2 T?p-n 



where a, stands for 



k,N 



2 



167r2 



as 



p—n 



+ higher twist 



and the An are the spin-averaged matrix elements of the spin-N twist-2 operator 



(5) 



(6) 



(7) 



The inclusion of quark charges is implicitly understood. The flavors of the quark-operators ip 
are combined to yield the proton minus neutron matrix element. {/J.i'} indicates symmetric and 
traceless combinations. The higher twist corrections are given by matrix elements of twist-4 
operators and were derived for the second moments of F^ and F2 in ref. 0. 



dx 



p—n 



4Q2 



p—n 



8g2 



p—n 



(8) 



3 



where we used the conventions of 0]. In this equation we only retained the twist-4 corrections. 
The target mass corrections are not exphcitly shown. Due to the high dimension of the 
operators it is at present not possible to perform such a calculation reliably in the framework 
of lattice QCD or QCD sum rules. An additional problem of such a calculation is that a 
renormalization scheme has to be found in which quadratic divergences in twist-4 matrix elements 
do not produce mixing with lower dimension twist-2 operators. That is still an unsolved problem 
in lattice calculations 



2J] . On the other hand, state of the art calculations of higher twist- 
corrections never claim an accuracy better than 30-50%. Therefore we claim that calculating the 
renormalon ambiguity in the coefficient function CL_Ar(Q^//i^, as) instead of the true higher-twist 
corrections to the longitudinal structure function is an legitimate procedure. The advantage is 
that such a calculation can be done for all therefore allowing to estimate the twist-4 corrections 
as a function of Bjorken-x. Note that a renormalon ambiguity in the coefficient function of the 
twist-2 spin-2 operator will only account for twist-4, spin-2, twist-6, spin-2 etc. operators and 
not for power suppressed twist-2 operators. This implies that target mass effects can not be 
traced by IR-renormalons. 

The truncated perturbative expansion of the coefficient functions of the moments of Fl and 
F2 can be written as 



mo(N) 
n=0 



(2) 

k,N Q4 



-2C 



(9) 



where we have accounted for the asymptotic behaviour of the perturbation series which makes 
only sense up to a maximal order mo(A^) depending on the magnitude of the expansion parameter 



With the standard normalization one finds B 



(0) 
L,N 



0, 



1 and B 



(1) 

L,N 



4Cf/{1 + N) 



Cf = 4/3 is the eigenvalue of the Casimir operator of the SU{3) colour group in the 
fundamental representation. We have indicated the ambiguity of the asymptotic expansion by 
including power suppressed terms. A^e~*" is a renormalization scheme independent quantity. 
C = —5/3 corresponds to Aj^g. We will show that only ambiguities up to order will 
appear in the 1/Nf approximation for the longitudinal coefficient function. Since non-singlet Fl 
and F2 are to leading twist accuracy determined by the same operators we obtain from Eq. (|^) 
with = ij?, as = as{Q^) 



mo{N) A2„-C 

V a^ 4- r(') ^ 



n=l 



mo(N) 



1+ E B^%a- + C^l^'-'~" 



n=l 



Expanding the denominator we find 



B^^^ a 



+ (4% - B^^B^l) al + 0{a'^ 



M2,n{Q') 



(10) 



{c'i^N - + C^^'lv) as + 0{al))]M2AQ') (H) 



Since the Callan-Cross relation gives -B|°7v — ^ ^^e calculation of the ambiguity in the 

perturbative expansion of Cl^n alone is sufficient to determine power-suppressed contributions 
to Fl~"'{x,Q^) up to 0{as/Q^) accuracy. 
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To extract the renormalon contribution to Cl,n we calculate the coefficients to all orders in 
tts in the 1/iV/ expansion where Nf refers to the number of active flavours [ffl, 0, ITBl. Each 



coefficient -B^^ can be written as an expansion in Nf 



= B't^ + B^r^Nf + . . . + Bi^^NJ'-^ (12) 

where the coefficient B^^ is unambiguously determined by the diagrams with one gluon-line 
that contains m — 1 fermion bubble insertions. These diagrams can be calculated comparably 
easy while the non-leading terms are much harder to evaluate. The non-leading terms are 
approximated in the procedure of "Naive Nonabelianization" (NNA) [jT9[ where the highest 
power of Nf is substituted by iVj ^ iV/ - 33/2 = -|/3o. Here 

/3o = 11 - ^Nf (13) 

is the one-loop coefficient of the QCD /^-function. The exact coefficient is therefore approximated 
as 

In cases where the exact higher order results are known, NNA approximates the exact coefficients 
well in the MS'-scheme In what follows we are going to calculate the coefficients B^/"^ of 
the {f3oO's)"^~^ expansion. For that we split the exact coefficient into 



^L,N - ^L,N + Ol,N ■ ll5j 

While B)[^^ contains only the effects of one-loop running of the coupling to order m, only S][ 
requires a true m-loop calculation. It will be checked a posteriori by comparison with those 



(m+l) _ 5(m) j-(m) 
L,N — ^L,N "T 'JL,N 

eric 

coefficient that are known exactly whether the neglection of 6^^^ is justified (see Eq. (pT]) and 
Table (|I])). Note that B^^^ = eJ^Jv = 4:Cf/{N + 1). In the following the NNA approximation 
to the coefficient function Cl^n{(^s) will be written as CL,N{cis). 

A convenient way to calculate Cl^n{0's) is to deal with its Borel transform 

The advantage of that representation is manifold. The Borel transform can be used as generating 
function for the fixed order coefficients 



(17) 

s=0 



and the sum of all diagrams can be defined by the integral representation 



1 °° 

- [dse-'/^^-'^^^BT[CL.N{asm. 
1^0 i 



Cl,n{cIs] „ 

Po 

Technically the most important point is the simplification of the calculation of the Borel trans- 
form of diagrams with only one fermion bubble chain. In that case the Borel transform can be 
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applied directly to the effective gluon propagator which resums the fermion bubble chain. The 
effective (Borel-transformed) gluon propagator is 



5T[a,D^^(A;)](s 



(19) 



In fact one only has to calculate the leading-order diagram with the usual gluon propagator 
substituted by the above one in which only the usual denominator of (— A;^)~^ is changed to 

To obtain the coefficient function Cl^n{^,cis) we have to calculate the 0{as) correction to 
the Compton forward scattering amplitude with the effective propagator Eq. (p^Q]). We get 



Cf 



T{s + N) 



(2 - -s)il + s + N) r(l + s)T{N) 

o r(s + Ar) 



r(i + s)r(iv) 



+ 



{2 + N){l-s) {3 + N){2-s) {2 + N){3 + N){1 + s + N) 



(20) 



The Borel transform exhibits IR-renormalons at s = 1 and s = 2. The position of the UV- 
renormalon s = —1 — N depends on the moment N one is dealing with. Formula Eq. ( 120|) has 
been derived independently in [p7| . 

The NNA approximants to the coefficient function in all orders in can be derived setting 
/i^ = and C = —5/3 in equation Eq. (pO|). It is interesting to compare the approximation 
of the NNA procedure with the exact results derived by Larin et al. for the non singlet 
moments N = 2, 4, 6, 8, denoted by Cl^^- 



CL,2(l,as 



C 



L,2 



C 



LA 



I, as 



c 



Lfi 



Cr 



a, ■ 1.77778 + 0^(74.963 - 4.54321Ar^) + a^(3238.62 - 392.56A^/ + 11.8957ArJ) 
a, ■ 1.77778 + 0^(56.755 - 4.54321A^j) + a^(2544.60 - 421.69Ar^ + 11.8957A^| 



N 



f 



-23.21 Y: qj 

as ■ 1.06667 + 0^(56.32 - 3.41333A^/) + 0^(2964.52 - 359.336A^/ + 10.889iV 



as ■ 1.06667 + af (47.99 - 3.41333A^/) + < (2523.74 - 383.052iVj + 10.889A^; 

Nf 

-15.18^9/) 

/=i 

as ■ 0.761905 + a,(44.4789 - 2.69569Ar^) + a^(2578.8 - 312.582Ar^ + 9.47219A^J) 
as ■ 0.761905 + a,(40.9962 - 2.69569A^/) + (^2368.2 - 340.069A^/ + 9.47219A^^ 

Nf 

-11-12 E?/) 

/=1 

2.23147A^/) + a^(2269.79 - 275.126iV/ + 8.33715ArJ) 
2.23147A^f) + a?f2215.21 - 305.473A^f + 8.33715A^" 



a, -0.592593 + 0^(36.8193 
as - 0.592593 + a^(35.8766 
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Exact results 28] 


NNA approximants 


N = 2 


Nf = 3 
Nf = A 


43.1254 a] + 1386.59 
38.5822 al + 1032.7 


61.3333 al + 2168. 
56.7901 al + 1858.71 


N = A 


Nf = 3 
Nf = A 


37.75 al + 1472.58 
34.3367 al + 1155.64 


46.08 al + 1984.51 
42.6667 0^ + 1701.4 a^ 


AT" n 

N = b 


Nf = 3 
Nf = 4: 


32.9091 + 1433.24 
30.2134 + 1152.07 


36.3918 a^ + 1726.31 
33.6961 al + 1480.03 


N = 8 


Nf = 3 
Nf = A 


29.1822 + 1373.83 al 
26.9507 + 1117.97 


30.1249 0^ + 1519.45 al 
27.8934 + 1302.68 al 



Table 1: Comparison of the NNA approximants to the exact results obtained in [|28[] for the coefficient 
function jv(l,as) up to order 0{al). We have omitted the 0{as) corrections which agree exactly. 



Nf 



•74 E?/ 

/=i 



(21) 



The sum over the quark charges Yl!fLi Qf stems from the exact calculation of the so called 
light-by-light diagrams where the photon vertices are connected with different fermion lines. 
Those diagrams first appear at three-loop level. 

The subleading Nf coefficients approximate those of the exact expression in sign and magni- 
tude. The leading Nfal and Njal coefficients of course agree exactly. The numerically important 
cases Nf = 3 and Nf = 4 are given in table 0. For comparison we have also given the NNA 
approximants to the exact 0{as) corrections to the coefficient function of the structure function 
F2 that were calculated in p9[. It is interesting to observe that NNA approximates the higher 
moments consistently better than the lower ones and gives better results for Fl than for F2. 
These features can be understood as follows. The most problematic property of NNA is the ne- 
glect of multiple gluon emission. As such processes are important for small x we cannot expect 
our NNA structure functions to be correct in this region. Ever higher moments of the structure 
functions are less and less sensitive to their small-x behaviour and therefore the NNA should 
systematically improve. 

As can be seen from Eq. (^) the perturbative expansion of Cl,n is not Borel summable. The 
poles in the Borel representation at s = 1 and s = 2 destroy a reconstruction of the summed 
series via Eq. (p^S]). Asymptotically the first IR-Renormalon, i.e. the pole at s = 1 will dominate 
the perturbative expansion giving rise to a factorial growth of the coefficient 



lim B 



(m) 
L.N 



am 

~ Po 



ds^^l- s 



(22) 



s=0 



This means that a perturbative expansion at best can be regarded as an asymptotic expan- 
sion and the expansion makes sense only up to a maximal value m = itlq ~ log((5^/A^). For 
higher values of m > uiq the fixed order contributions will increase and finally diverge. The 
general uncertainty in the perturbative prediction is then of the order of the minimal term in 
the expansion. It can be estimated either directly or by taking the imaginary part Q'/vr (divided 
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Exact results 29 


NNA approximants 


= 2 


Nf = 3 

Nf=4: 


1.69381 al 
-3.63952 al 


71.9999 
66.6666 


N = A 


Nf = S 

Nf=4: 


91.3793 al 
74.3914 al 


229.3403 
212.3488 al 


N = 6 


Nf = 3 

Nf=4: 


218.3608 al 
190.3477 a'i 


378.1762 
350.1631 a2 


N = 8 


Nf = 3 

Nf=4: 


357.0326 
319.1078 a'i 


511.9848 al 
474.06 a'i 


N=10 


Nf = 3 
Nf = 4: 


498.6271 a'i 
451.7658 a2 


632.6276 
585.7663 



Table 2: Same as Table |I| for the coefficient function C2,Ar(l, cis) of the moments of the structure 
function Ff^". These were obtained in up to order 0{aj,). On the right column we compare 
these with the NNA approximants. The 0(a<j) corrections agree exactly. 



by vr) of the Borel transform [|TT]. From Eq. (|T8|) we get for the function Eq. (|2 
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^Cf Ace-^ N ^ ACf A^e"^^ A^^ + 



f3o N + 2 (3o N + 3 



(23) 



The ambiguity in the sign of the IR-renormalon contributions is due to the two possible contour 
deformations above or below the pole at s = 1 and s = 2. For the moments of we then get 



Observing that 



1 + N 



\ A2_e5/3 



(24) 



dx x^-^ [6{x - 1) - 2x^) 



N 







A^ + 2 



the above equation is easily transformed from the moment-space to Bjorken-a; space. 

Fnx,Q') + ^^Fr''ix,Q')=ACFa.iQ') j ^ Fr"(y, Q^) + 0(a^ 



^CF^rsfl 



Fr"(x, g^) - 2 / ^ (-] ' Fr"(y, Q') + Oia, 

■I y \yJ 



Q i y 



(25) 



(26) 
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We have neglected the contribution of the second IR-renormalon since it is of the order of 1/ (Q^)^ 
while there is a contribution of order a^/Q^ related to the ambiguity in the coefficient function 
of F2 which we have not included. We have included the kinematical and target mass corrections 
to the order we are working as given in 

In connection to the IR renormalon, it is interesting to investigate the corresponding ambigu- 
ity in the definition of the twist-4 matrix elements. This can be done for the second moment of 
Fl where the contributing twist-4 operators are known, see Eq. (P). Composite operators have to 
be renormalized individually and have their own renormalization scale dependence. Operators of 
a higher twist, and therefore of a higher dimension, exhibit power-like UV divergences in addition 
to the usual logarithmic ones. In particular, a quadratic divergence of a twist-4 matrix element 
contributing to Fi, results in a mixing with the lower- dimension twist-2 matrix element. When 
the calculation of the matrix element is done in the framework of dimensional regularization, the 
quadratic divergence does not appear explicitely. It manifests itself as r(l —d/2) factor, singular 
at (i = 4, and at c? = 2, where it corresponds to usual logarithmic divergence. Evaluating the one 
loop contribution to the quark matrix elements of twist-4 operators, with the Borel transformed 
propagator (Eq. (|l^)), in 0? = 4 dimensions we obtain an expression singular at s = 1. 



BT 

BT ^ 



C^^iasMi' (D'^gG.^,) 7.}V|P)"-"] (^)|_^ = ^ f ^ ) {P\H,^D,}^\P) 



\p—n 



(27) 



1 - s V 
28 Cf / /i"e 
Yl-s \Q 



The singularities at s = 1 is the manifestation of quadratic UV divergences. Inserting this 
together with Eq. ( [20| ) into the operator product expansion for the second moment of Fl we see 
that this contribution indeed cancels against the IR-renormalon contribution to the coefficient 
function of the twist-2 term, resulting in an expression which is free of perturbatively generated 
ambiguities up to the A^/Q^ order. 

The appearance of an IR-renormalon ambiguity in a perturbative calculation thus indicates 
the need to include higher twist corrections to interpret the perturbative expansion to all orders. 
Of course one can argue that the numerical value of the IR-renormalon uncertainty has no 
physical significance since it has to cancel in a complete calculation. As we explained above, in 
some case a higher twist estimate based on IR renormalon ambiguity has proven to be a fairly 
good guess, at least for low moments of the structure functions. In the present case the IR 
calculation can be easily done for all moments, so that the result can be extended to produce 
a model of the the full x-dependence of higher twist effects. We keep in mind that such a 
model cannot have significance beyond phenomenological level for the following reason. The 
renormalons ambiguity in the coefficient function Cl,n{Q^) is a target independent quantity 
of pure perturbative nature, while "genuine" higher twist matrix elements are a measure of 
multiparticle correlations in the target and are process dependent. 

It is interesting to compare our prediction for the twist-4 part of Fl with the available 
experimental data. To this end we use the phenomenological parametrizations for R{x, Q"^) 
Eq. (^ of 0. To extract Fl{x,Q'^) we have choosen the parametrization of F^a:, Q^) and 
F^{x,Q^) of Q valid in the region 0.006 < x < 0.9 and 0.5 < < 75 GeV^. With ~ 



2(^1 - Fi) we have 



FL + -^F, = ^^^ + -^\F, = Ff-^-\x,Q') + ^^ + 0\-]^ (28) 
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Figure 1: Comparison of d'^'^^{x,Q'^) (long dashed line), the power suppressed contribution on the 
right hand side of Eq. (P^, with the phenomenological fit d^^{x,Q^) (full line) Eq. (|2B|). We have 
also plotted the IR-renormalon (twist-4) part (short dashed line) and the target mass corrections 
(dotted line) separately. The agreement with experiment is much better including twist-4 corrections 
than without them (dotted line). We have chosen Ajjg = 250 MeV, = 5 GeV and Nf = 4. 

where we neglect the 0{1/Q^) contributions for consistent comparison with our calculation. 
The dependence of the higher twist coefficient Q^) is only logarithmic. In Figure 1 we 

compare the experimental fit of the higher-twist coefficient with the Q CD-calculation (Eq. (p^)), 
where we have shown target mass and twist-4 contributions separately. We observe a rather 
large contribution coming from the IR renormalon estimate for the twist-4 part, which accounts 
for more than half of the discrepancy between the experimental fit and a prediction which takes 
into account the target mass correction only. The sign of the IR renormalon contribution, which 
cannot be determined theoretically, should be choosen positive. This leads to an astonishing 
agreement (possible somewhat fortuitous) with the experimental fit. Our estimate based on the 
calculation of the IR ambiguity has proven to be phenomenologically surprisingly successful, 
predicting a high twist-4 contribution to Fl in accordance with experimental results. It further 
supports the idea that, while the rigorous QCD calculations of higher twist contribution to 
are not yet available, calculations like the one presented in this paper can be used to predict the 
order of magnitude of power suppressed corrections. 
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